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ABSTRAGI 


An analysis is made of current Navy procedures for 
assigning a probability density function to demand and 
the technique(s) used to forecast the parameter(s) of 
the particular density function chosen. The gamma den— 
sity function is investigated as a possible replacement 
for the three density functions currently being used. 

A comparison of the gamma and normal density functions 

is made with regards to inventory costs, observed protec 
tion levels, and unit effectiveness. Additionally, a 
comparison is made, through a simple simulation model, 

of demand patterns generated by demand reporting as it 

is done today and consumption data reporting as it might 
be done in the future. The intent of the simulation is 
to provide some insight into the impact the demand re- 


porting method might have on the variance of demand. 
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I. INTRODUCTION 
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Today, the Naval Supply System faces the task of sup—- 
porting very complex weapon systems under the burden of a 
funding climate that is already austere and will probably 
become more so. Accordingly, this situation makes effi- 
cient inventory management of utmost importance. This is 
not to say that inventory management has not been vital in 
the past, but rather the decision-making process with re- 
gards to inventory management becomes more critical. Ob- 
viously as the dollars the supply system receives are 
reduced, the application of the dollars must be judicious— 
ly applied to provide the best support possible to the 
fleet. By support it is meant the repair parts required 
to keep a weapon system or essential shipboard/aircraft 
operating. In view of the fact that an effective essen- 
tiality coding has not been established within the Navy, 
this means supporting almost everything with a Navy stock 
MUD Ole, Shout be noted that the general objective of 
the Naval Supply System Command for the past several years 
has been to reduce the time equipment is down or not 
operating due to the failure of a repair part. Theses 
course entails more than having the part in inventory; 
however, this is the most critical aspect since most items 
required are unique to the Navy and are not readily avail-— 


able from industry. 





Essential to having the right part On hand at the right 
time is the ability to predict future demand for repair 
parts. This 1S by no means an easy task in view of the 
random nature of demand over a given time interval. The 
purpose of this thesis is to appraise what is being done 
within the Naval Supply System to predict future demand 
and to investigate possible alternative techniques. 

In predicting future demand there are two major factors 
that must be addressed. They are: 

1. Determine the probability distribution that repre- 
sents or at least approximates the actual distri- 
bution. 

2. Develop a technigue for forecasting the parameters 
of the probability distribution determined to be 
representative of how future demand is expected to 
be generated. 

The current Navy methods of forecasting demand will be 
investigated. The selection of the probability distribu-— 
tion will then be discussed with primary emphasis on com— 
paring the normal distributions presently in use with gamma 
distributions as to their effect on the establishment of 
inventory levels. A.brief discussion is made with regards 
Oo fhe use of the Poisson, negative binomial and the pseudo— 
normal distributions. The method of demand collection is 
then examined, and some possible problems associated with 
the current system within the Navy are discussed. Finally, 


the conclusions and recommendations generated by this study 


at CmwEc Seneca. 





ITI. DEMAND FORECASTING 


There are many ways to forecast demand; however, most 
of them can be grouped into two broad categories. The 
technical or subjective method and the scientific (some 
people would prefer mathematical) or the objective method. 
The technical method is the method by which the technical 
experts in a particular field give their best subjective 
estimate as to the demand for an item. This method appears 
appropriate only when there does not exist past demand 
data for an item. It could be useful possibly as a weight— 
ing factor in one of the scientific methods. Presently, 
the Navy supply system uses the technical method for com— 
puting Best Replacement Factors (BRF) for new items of 
supply or for items which have not realized any demand 
but are required to be stocked. This method will not be 
discussed any further, since the objective of this thesis 
is to deal with the items having recurring demand. 

The scientific method of forecasting demand breaks 
down to baSically two types of technique. The first in- 
eludes themtimesscricss techniquesmwheres past demand pat— 
terns are used to predict future patterns. Some examples 
of this type are the straight average, moving average, 
single and double exponential smoothing, spectral analysis 
and BayeSian procedures. These techniques rely solely on 


looking at the pattern of past demand over time to predict 








future demand for a specified time frame. The second type 
within the scientific method is the causal type. Although 
this type also relies on past demand patterns, it attempts 
to explain the causes for the particular patterns. MThere—- 
fore, this technique requires the knowledge of certain 
causal variables which in turn can be related to the depen-— 
dent variable for prediction purposes. 

This is the ideal technique; however, it 1s very time 
consuming, costly, and usually very difficult to determine 
the causal variables. This technique is usually related 
to econometric models. The Navy uses a regression model 
in determining the Mean Absolute Deviation (MAD) of demand, 
which is further used to determine the standard deviation 
of demand. In this model the explanatory variable or 
causal variable is the average demand D. This model will 
be discussed in more detail later in this section. 

The following time series techniques are defined: 

A. STRAIGHT AVERAGE — This technigue does nothing 
more than divide the total number of demands by 
the number of obServations as dictated by the 
desired time interval. For example, if average 
monthly demand was desired, the total number of 
gemands) wou ld® bese idecdaby ene: total number of 
months for which demand data was available. 

B. MOVING AVERAGE — In laymen's terms this method 
drops the oldest observation of demand and adds 
the latest and then computes the straight average. 
The computational form selects a number of base 
observations, say four months, eight months, or 
four quarters. Once the basSe is selected, a 
forecast is generated by computing the average 


of the base observations. Then when a new 
Observation 1S available, the new forecast 1S 





computed by adding to the old forecast the new 
observation minus the oldest observation divided 
by the number of base observations. Mathemati- 
Seplglkyy 


NEWEST OBS—OLDEST OBS 


NEW FORECAST = OLD FORECAST + Number of base obs. 


C. EXPONENTIAL SMOOTHING — This technique is used 
to assign different weights to different observa— 
tions of demand. In the straight and moving 
average methods equal weight was assigned to each 
and every observation. In exponential smoothing 
this is not the case. The more recent demands 
are given heavier weights than the older demands. 
Mathematically the forecast is calculated as: 

NEW FORECAST = OLD FORECAST + a(NEW OBS — OLD FORECAST) 
which can be further simplified to 

NEW FORECAST =(1—-a)OLD FORECAST+ (a)NEW OBSERVATION 


where qa = aSSigned weight to most recent observa— 
ede ye tebe 0) cs eh eI 


Now that the forecasting techniques have been defined 
it would be worthwhile to look at the advantages and dis— 
advantages of each technique. An advantage of the straight 
average is the fact that it provides the maximum likelihood 
estimate of the mean demand when a normal distribution de— 
scribes lead-time demand. In other words, the sample mean 
1S an unbiased efficient estimator of the actual mean of 
the distribution. The disadvantage of using the straight 
average 1S its slow reaction to large increases or decreases 
in demand. To use this latter point as an argument against 
Phesuse Crethnis technigqgueyseems te be Contradictory f€o the 
assumption of lead-time demand being normal with a constant 


mean. Since this technique produces the best estimator of 
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the parameter, it would seem that the more observations 
available the closer this estimate would be to the actual 
parameter regardless of the immediate consequences of a 
Significant change in a particular demand observation. 
Another disadvantage of the straight average lies in its 
inability to adjust to changes in demand patterns caused 
by a varying mean. 

The advantage of the moving average over the straight 
average is its ability to react more quickly to changes 
in demand patterns. The disadvantages associated with 
using this statistic are high computer costs in both space 
and time, which are due to the requirements to store all 
past observations over the base period. 

Exponential smoothing provides an estimate of the mean 
which is asymptotically unbiased and it reacts more quickly 
to changes in demand patterns. It also requires very 
little computer space and time. These latter attributes 
make it the most desirable technigue as far as the Navy 
Supply System is concerned. 

Presently the Navy Uniform Inventory Control Proce- 
dures (UICP) uses exponential smoothing for forecasting 
demand; however, this forecast is refined by use of trend 
tests and a tracking signal filter. The trend test is 
used to verify if there is actually a change in demand 
pattern. The tracking signal iS a statistic which is used 


to determine if the forecasting rule is performing well. 
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It is an algebraic sum of forecast errors. If the fore—- 
casting rule is operating properly the sum of thesSe errors 
should be close to zero; if it is large, then the system 
is considered out of control and recomputation is performed. 
So far, only techniques for forecasting the mean 
(average demand) have been discussed; however, the variance 
which iS a measure of how demand varies about its average 
must also be forecasted. Estimating the variance accurate— 
ly is critical since it has a pronounced affect on the 
reorder point and safety levels which of course drive the 
stock levels maintained. The stock levels further determine 
the major portion of the inventory investment cost. Pres— 
ently the UICP model of the Navy estimates the variance by 


using the mean absolute deviation (MAD). The MAD is defined 





as the expected value of X— | and the most commonly used 
estimate of MAD is the sample average of the absolute differ-— 
ences between the observed demands and their sample mean, 
|x.—x | 


N 





Lota = 


i1=1 


ire ULrCe modeledoes noe use this Formulation but rather 


B 


uses the formula, MAD = cp , where D represents expected 
lead-time demand. The formula is based on a study done by 
Fleet Material Support Office (FMSO) where regression analy— 
Sils wassused to determine the coefficients a and Bg. This 


study {l] was done in 1962, and the values of a and 8 were 


determined to be 1.37 and .717, respectively. 


2 








The UICF model takes the MAD computed as above and 
then uses as the estimate of the standard deviation 
Shr _ 25 MAD. Once a MAD is determined for an item 
SuccesSive MADsS are exponentially smoothed the same as 
the estimate of mean demand. The use of this technique 
for estimating the standard deviation is based on work 
done by R. G. Brown [2]. This technique has come under 
criticism by several investigators including particularly 
Dr. P. W. Zehna [3] when used with a probability density 
function other than the normal. 

Using a 1000 item sample of demand data provided by 
the Air Force, which has been used in studies done by all 
services in the past several years, a 50 item subsample 
was selected randomly and estimates of the standard devia- 


tion and MAD (MAD, = A baad form) were computed. The 


1 
sample standard deviation was then compared with the value 
oe (MAD, ) and a sign test performed to test the hypo- 
thesis that there was no difference between the two esti- 
mates of standard deviation. Since the difference between 
the two standard deviations was always positive, the sign 
test strongly rejected the hypothesis. 

sel view ©f£@ehe faectethatwenes translation between tlic 
theoretical MAD and the standard deviation appears to be 
questionable it seems prudent to also ive Stelga ae SEhesucc 


of the linear regression technique of estimating the MAD. 


Problems encountered with this technique would certainly 
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compound any problem associated with utilizing the MAD to 
estimate the standard deviation as is presently being done. 
Great care must be taken when uSing a linear regression 
prediction scheme where the values of the explanatory 
variable (the predictions of mean monthly demand) used in 
estimating the least squares coefficient are not known. 
Using a regression model to predict beyond the range of 


values of the explanatory variable can lead to serious dif- 


ficulties. From the random sample of 50 items, the regres— 
Sion eStimate of the mean absolute deviation, MAD = 1.37 
x fi? was calculated for each item. These values were 


compared with the other estimate of the mean absolute devia-— 
tion, MAD, . It was noticed that MAD was always greater than 
MAD, when MAD, was less than one and MAD tended to be less 
led ekcbel MAD, when the latter value exceeded one. The hypothesis 
that the two estimates of mean absolute deviation are the 
Same was tested uSing the Wilcoxon Signed Rank Test, and 
1t was strongly rejected at the 0.01 level of significance. 
Results of that test are found in Appendix A. 

The above results combined with the earlier comparison 
of the sample standard deviation and 1.25 MAD, meint “Out 
the possible source of problems experienced by the Navy in 
estimating variance. The wide disparity between the differ- 
ent estimates of the standard deviation can best be illus— 


trated by presenting the three estimates for each of the 


SUM EECnS was 1S cones i nm rable I, 
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TABLE! 
Comparison of three methods for estimating 


standard deviation of demand. 
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The method of forecasting the two critical parameters, 
the mean lead-time demand and the standard deviation of 
lead-time demand, discussed in this section will also be 
compared in Section III in order to appraise the effects 


on inventory leveis caused by their use. 
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Pi Sob LE@TION @F PROBABMAIRY DENSI 
FUNCTM@NS TO" REPRESEN BEAD—TIME 
DEMAND 
Whereas there has been a great deal of concern in the 
past several years regarding the ability to forecast lead—- 
time demand and variance, little or no concern has been 
evidenced regarding the demand Aree ton itself. Ob- 
viously no one probability distribution function is going 
to fit the demand for all military items; however, one is 
interested in determining if there is any one distribution 
that can best represent demand as the real world sees it. 
An affirmative answer to this question would require a ver- 
Satile distribution function since it is readily recognized 
that demand for military items covers the range from very 
low to very high, and frequently is erratic in nature. The 
second guestion to ask, given the first is answered affir- 
matively, would be the feasibility of using such a distri- 
bution. Three Saopaloi ey distributions are used in UICP 
in varying degrees we describe demand. These distributions 
are the Poisson for very low demand items (i.e. average 
annual demand < 1), the negative binomial for items with 
annual demand in the range of two to ten, and normal for 
items with average annual demand greater than ten. 
Prichard and Eagle [4] state that empirical demand dis— 


tributions for items with low means are usually skewed to 


the right thereby making is desirable to use the Poisson 
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and negative binomial distributions to describe this type 
of item. They further state that the Poisson, negative 
binomial, and normal distributions seem to satisfactorily 
approximate demand for a majority of items and also indi- 
cate that price and consumer behavior usually are the 
causes for differentiating between the use of the Poisson 
and the negative binomial. That is, a cheap low demand 
item will be ordered in batches thereby cauSing a high 
variance and suggesting the use of the negative binomial. 
On the other hand, an expensive low demand item will be 
ordered one at a time thereby being better represented by 
the Poisson. Although not specifically stated, the reader 
1s led to the assumption that the high demand items can be 
represented by the normal regardless of the price. It is 
by this reasoning that they contend a single mathematical 
distribution cannot represent all demand. They do not 
address, however, the possibility of a single distribution 
Satisfactorily approximating these three distributions. 

As mentioned previously, the Navy Supply System is using 
these distributions in varying degrees. A recent study 
completed by FMSO indicates that one Navy Inventory Control 
Point (ICP) is using the normal to represent all of its 
demand [{5]. Additionally, the use of the Poisson has de- 
creased in the past several years, as tables have been 
generated to use in its place. Basically when discussing 
demand distributions within the Navy it is the normal 


distribution that is being referred to. 
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Recently there has been a move underway to substitute 
the pseudo—normal distribution for the normal and the 
other distributions. The pseudo-normal has been recom— 
mended because of the computational advantages that it 
allows in the calculation of inventory reorder levels and 
reorder points. Because it performs about the same as 
the normal with regard to describing demand patterns, we 
will concentrate our attention on the normal distribution 
in this study. 

As stated previously, representing military demand by 
a Single probability distribution would require a very 
versatile distribution. It is felt that the gamma family 
of distributions best meets this qualification. Dr. 
Peter"W. Zehna [6] states: "The family (gamma) 1S so ex— 
tensive in Shapes of densities available that it isa 
fairly safe assumption to make aS a model for an experi- 
ment described by almost any non-negative random variable." 
FOGdaastonc bamnOc Lymer Ort, sandmsotone (gi point Out that an 
most cases involving a random variable X which is known 
to be positive the assumption that X has a gamma density 
will provide an approximation or at least an inSight into 
the true but unknown Situation. Finally, the versatility 
of the gamma is shown graphically in Figure 1 where curves 
associated with the Poisson, negative binomial, and normal 
distributions are displayed in 1A and the curves in 1B 


represent a particular set of gamma distributions obtained 
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by varying the parameters A and a. These parameters re- 
ferred to are called the shape (A) and scale (a) para- 
meters of the gamma function and will be discussed in 


more acta ial later. 


£ (x) £ (x) 


=. 


A x B Did 
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A true but unknown situation seems to describe more 
than adequately the situation inventory managers find 
themselves in when trying to predict how demand is going 
to occur. Therefore it seems appropriate that the gamma 
PHowha” Be@inuvestigaecasds EQ tseability to represent 
demand of military items and to be compared with the normal 
ais tO which better diescribes the situation. It 1s pre- 
cisely this objective that the remainder of this section 
1s devoted to. 

Using the sample of 50 items chosen from the 1000 item 
sample of Air Force data, the monthly mean Ce Standard 


deviation (S) and mean absolute deviation (MAD) were 
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computed based on 57 months of past demand. Additionally, 
the gamma parameters A and a were computed using the 


method of moments in the following manner: 


FOR I[(x,A,a)= aay Ox ee E(x)= A/a, V(x)= A/on 


THEN solving for A and a 


x/s7 = a and Sisk = A 


These parameters were computed mechanically by the computer 
which at the same time was performing a goodness of fit 
test for two different cases. In the first case the normal 
distribution was fit to the data using the sample mean and 
the sample standard deviation which are the maximum likeli- 
hood estimators for the parameters of the normal distribu- 
tion. In the second case the gamma distribution was fit 

to the data using the sample mean and variance to compute 
the estimates of the gamma parameters A and a. A Kolmogorov— 
Smirnov (K-S) test was performed in each case resulting in 
a strong rejection of the normal hypothesis and an accept— 
ance of the majority of items under the gamma hypothesis. 


Table II provides a summary of the results of the K-S test. 


TABLE II 


Results of Kolmogorov—Smirnov test on 
50 randomly selected items. 


Distribution No. Items Ttems Percent of Items 
Fitted Accepted an ccin Accepted 


Nomamees | 0 50 03 
Gamma (A,qa) oF 13 74% 
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A listing of results obtained in both cases is provided 
in Appendix B. 

The necessity of fitting a probability distribution 
SemeCniakdsdcpameone SmimOm stbe Lequirementweto protect 
against running short of stock. In inventory models pro- 
tecting against shortages is achieved by ordering stock 
before the amount on hand falls too low. The procurement 
policies specify that orders will be placed as soon as 
the stock assets reach or fall below a given level called 
the reorder level. The higher the reorder level, r, the 
lower will be the probability that the total demand dur- 
ing the lead-time exceeds r. Define the risk to be this 
probability and let its complement be called the protection 


level. 
RISK = Pr [Total lead-time demand > r] 


The reorder level is adjusted up or down with the degree 
of adjustment depending on the desired protection level 
and the probability distribution of lead-time demand. fThe 
amount of ieeRe Reel ordered, the reorder Gvdantley, «ss usually 
constrained by the administrative cost to place an order, 
EDS Cosi ee nolans thai aAnOulnemor mitt eCilawe lm LnVentLony 
and the cost to’the system of being out of stock when a 
CON tC ocn th Se OeCSSecOSESeaLe Notmally referred to as 
ordering, holding, and shortage costs and are not easily 
obtainable. Consequently they are usually a product of 


the particular inventory system. 
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The Navy uses a modification of Hadley and Whitin's lot 
size-—reorder point model (<Q,r> model) defined as follows: 


i LGEsfsbe Lime annuals holding cost rate per unit of 
stock, Ax the fixed cost of placing an order, C 
HiceecOsE pCOsuMmtE fOr 1bem ie te mean annual demand 
FOG item 1,.and 7 thegshortage cost for each unit 
short. (A fixed shortage cost is not assigned for 
each item; rather, the penalty 7 is simply mani- 
pulated to meet budget and risk constraints.) If 
H(r) 1S the probability that the total lead-time 
demand for item 1 will exceed r, the reorder level 
is determined from the equation 


LC. 0 
C5 ns) =, Se 
es Ose Tt 
1”i 1 
ee PmC=reOnaer GWanltity, Ol," fer item 1 is taken _to 
be the maximum value among (1, Qi, p./*) - Where 


‘i 
QO. _ Se 
kee 

Lk 


is the well known Wilson lot size. 

It 1S important to note that under the Navy UICP model the 
computation of Q. is completely independent of the reorder 
point and the risk, thereby making Q. independent of the 
demand probability distribution. 

ae Pemc Ola Hatten Bo MPCP enGdc Reinder =the aSSump= 
tions of normal versus gamma distribution of lead-time 
demand it iS necessary to obtain a feeling for the rela-— 
tive difference between the theoretical reorder points 
assigned by the gamma versus those assigned by the normal. 
For an illustrative purpose three items were chosen from 


the 50 item sample. Table III shows the results of 
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comparing the empirical, gamma, and normal distributions. 
Appendix D gives the graphic representations of the 
empirical cumulative distributions from which the actual 
protections were determined. The two normal distribu- 
flons dmtter onlysim the choice of standard deviations. 
In one case the sample standard deviation is used and in 
the other case SD is computed using the regression for- 
mula for the MAD as discussed in Section II. The two 
sets of figures for each distribution represent 80 and 50 
percent protection levels respectively. Table III indi- 
cates the gamma does a much better job at approximating 
the empirical distribution for the first two items than 
does either of the normal distributions. For the third 
item the gamma, and the normal using SD, both do a good 
job of approximating the empirical distribution at the 80 
percent protection level. At the 50 percent protection 
level only the normal with SD does well. At this level 
the gamma underprotects as much as the normal uSing S over- 
PrOtCeECts wee LOOking atlas hiistogram of this item (Figure #2) 
it can be seen that the gamma is not sensitive enough to 
the three large monthly demands (662, 757, 1733) and the 
normal 1S too senSitive. This situation seems to be the 
crux of the inventory problem. Should the system over 
react to peaks in demand or should it ignore them? Ob- 
viously it costs inventory dollars to over-react such as 


in this case of stocking 71 more units than required; 
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however, how does one measure the cost of the shortages 
one can anticipate by only stocking 16 units. It is in- 
teresting to note the effect the computed standard 
deviation, SD, has on this item. The value for SD is 
approximately 20 percent of the value S. This is not 
difficult to understand when one considers the fact that 
in the computation of the MAD the mean is raised to the 
-717 power and then multiplied by 1.37. Except for mean 
monthly demands less than 3.04 this formulation will al- 
wayS give a standard deviation less than the mean. 

The final comparison between the gamma and the normal 
is made with regards to the inventory levels and the 
associated costs that are generated under each distribu- 
tion. A FORTRAN program was written to compute reorder 
levels and reorder quantities for each of the 50 items 
chosen from the demand data base. These values were then 
used to generate holding costs. In addition, shortages 
were accumulated and an overall unit effectiveness was 
computed. Overall unit effectiveness is defined as the 
total number of units satisfied divided by the total num-— 
ber of units demanded. Appendix C contains three programs 
to accommodate the computation of the gamma cdf, and the 
use of the exponentially smoothed average and the predicted 
Standard deviation for the normal. Two different runs 
were made to determine the impact that different parameters 


would have on the levels computation and associated costs. 
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The procedures of the program for the first run are as 
follows: 


1) Read the following values: 
Holding Rate = XI = .01/month 
Ordering CostsFAx*=* $25. 


2) Generate unit price by randomly selecting values 
between one cent and 50 dollars. The same seed 
for the random number generator was used to 
obtain the same sequence of unit costs each time 
the program iS run. 


3) Compute the means, variances, and standard devia- 
tions using the first 21 months of demand history. 


4) Compute the initial reorder levels and quantities 
and randomly set the initial assets on hand between 
r+] and r+Q. It 1S in this step where the com— 
PUitduLomeo: eb itSattected by both the probability 
distribution and the parameters chosen and the Q 
By the particular forecasting technique chosen 
(1.e. sample mean, exponential smoothing, etc.). 


Pee eine ne PrLoglan tor tne Lehalnmimg so months Conpar= 
ing each month's actual demand with the available 
stock on hand and computing units snort, if any, 
total holding cost and finally overall unit effec— 
tiveness. 

6) The final step 1S to compute the observed protec— 
LIOhMe LaAlS 1S domes by divading the total number 
©f Jead—time periods adn whichishorts do not occur 
by the total number of lead-time periods. 

The first caSe was run three times for protection levels of 
SOs, 20 Percent f£OG Cach elLwehe Ehree™ distributions; 

the normal uSing the sample mean and sample standard devia- 
tion; the normal using the*exponentially smoothed average 
and the exponentially smoothed estimate of the standard 
deviation; and finally the gamma using the parameters estli— 


mated by the method of moments. The results are provided 


an Tables iv. 
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TABLE LY 


“OBSERVED PROTECTION, TOTAL HOLDING COST, TOTAL UNITS SHORT, 
AND OVERALL UNIT EFFECTIVENESS FOR VARIABLE r and Q = MAX 


(1.0, OQ, t/4). 


© 


% 


Des tred=Pretection bevc a—-s 


Normal Gamma Normal 
(x,S) (A, 0) (SMO,SD) 
Observed 
Protection (A se shell Wo. 67% er aoe 


Total Hold- 
ang Cost oe oe 2 cee ae wr 3) 


Rotea Lani ts 





SiOre i, om Lage gl 
Overall Unit 

Effectiveness Soe = e8./4¢ 

Desitrecomrrelectie@em Levcls——o0s 

Observed 

Protection Oak S ‘Sigilypggulgis > 56.90% 
Total Hold— 

ing Cost 7 OTS Ae oo 61, 236 


Total Units 

see Wir re 
Overall Unit 

Effectiveness! 93.44% 3 (RSE: 92. 14% 


Desired Protection Level = 


Observed 

Protection 59.388 50°. 46% 70.003 
Tosca HoLd— 

ingiliest amhe ery 
Total Units 

ai oe 409 
Overall Unre 

ahh oy ian egy ee oie 96.49% 
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Significant among the results depicted in Table IV are 
the high unit effectiveness corresponding to the low ob— 
served protection levels. In fact the observed protections 
fell short of the desired levels in all but two cases where 
the normal distribution was used with the parameter esti- 
mates SMO and SD at the 30 and 50 percent desired protec 
tion levels. These resultS point out that the system can 
satisfy most of the demands (high unit effectiveness) while 
achieving very low protection levels. Both the short lead— 
time (taken to be one month in this evaluation) and the 
large values of the reorder quantities Q relative to the 
reorder levels r account for this phenomenon. For in these 
circumstances the system rarely relies on the reorder level 
to protect against stockouts. In most cases the available 
stock iS greater than the reorder level. 

Both the protection level and the unit effectiveness 
figures obtained by the gamma distribution fall below the 
figures obtained by the normal distributions. At first 
this seems an indictment against use of gamma distribution. 
However, one must consider the extra cost of achieving this 
higher protection. That is, one must ask if the trade off 
1s cost effective. 

To obtain higher protection levels the normal distri- 
butions must be choosing higher reorder levels on the 
average. These higher protection levels result in higher 


investment and holding costs. Thus a comparison of the 
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holding costs gives an indication of the cost of achieving 
the higher protection. Although no definite conclusions 
can be drawn from the values in Table IV there are indica-— 
tions that the gamma performs better than the normal with 
the same estimates of the mean and standard deviation at 
the higher protection levels when measured in terms of unit 
effectiveness per dollar of holding cost. Also, when com— 
paring the two normal cases, the one with the exponentially 
smoothed estimate appears to perform better than the other 
uSing unit effectiveness per dollar of holding cost. 

In order to better assess the effects of the type of 
probability distribution on the measures of effectiveness 
and the costs, the reorder quantities were all forced to 
be Q=l1. This allows us to focus more clearly on the impact 
of the reorder levels determined by the two types of dis-— 
EY TOWe1On . “Also only the two cases where the sample mean 
and the sample standard deviation were used to estimate 
the parameters were considered so that any differences 
eed Ga bomecenes  DUteCdmEOe thewt _Eemenmene probabality drstri— 
butions. Table V shows the results obtained with the re- 
order quantities set at one. Only the obServed protection 
levels and the holding costs are presented in Table V 
because when Q=l1, the protection level should be a better 
measure of performance than unit effectiveness. 

The dominance of the gamma over the normal is better 


illustrated by Figure 3. There it 1s clearly shown that 
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TABLE V 


OBSERVED PROTECTION AND HOLDING COSTS 
Mie shLOUomoORVELO OF PROTECTION: 


Desired Protection Level 303% 
Demand Distribution Normal (x,S) Gamma (A,a) 
Observed Protection (OP) PAS pp Os Zoe S 
ToGolmHoOldaing Cost (THC) eck ies Jo 
503 
OP 62.53 33.6% 
THC J tev S\45 2265 
60% 
OP 68.73% 42.9% 
THC 19,868 SRG ay ; 
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OP 75.52 50.25 
THC 29,018 By 8) GS 
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the gamma protection is greater than the normal protection 
for any value of holding cost. These results coupled with 
the comparisons involving unit effectiveness versus dollars 
of holding cost support the hypothesis that the gamma is 
more cost effective than is the normal. 

The values obtained for actual protection in Tables IV 
and V appear to be inconsistent Since one would expect 
greater protection with higher inventory levels. The prob- 
lem appears to be in the definition of protection and 
relates to the frequency at which the system reaches its 
reorder point. The effect that Q can have on the protec- 
Eon level is best shown by an example. A year's demand 
for one of the items of the 50 item sample 1S recorded in 
Table VI. or Gase Of computation the variable Q was com— 
puted as the total yearly demand divided by 4(Q=t/4) and 
the reorder level was calculated from a normal distribution 
Pom vema 0 mercent protection level. The initial on hand 


level of stock was assumed to be rl. 


TABLE VI 


PERFORMANCE COMPARISON OF Q=1 vs. Q=t/4. 


Demands 4 Soe ie yee 45 Oe Lae 
Q=T/4=51 x= 17 
Ona way abe 52 38 —7/9*% =—20% 24 19 19 19; r= 17 

Protection = 1 — 2/3 = 333 OH=18 
U.E. = 1 — 99/203 = 55% 

QO = 1.0 

OH 14* 10*15*18 16* 4*—94* —25* 11*13%*18 18 

Protection = 1 — 2/9 = 77% 
U.E. = 1 — 124/203 = 39% 


where * implies an order was placed; U.E.=unit effectiveness 
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Finally, a discussion of the overall unit effective~ 
ness figure used in Table IV 1S warranted. An important 
question being asked today is whether an inventory system 
should be geared to protect against total units short or 
the total number of requisitions short. Obviously, if 
only one unit were demanded per reguisition (which is an 
assumption of the Q, r model) these two measures would be 
the same. However, when examining demand patterns for 
the items used in this study it was noticed that there 
would be many months of demand within a particular range 
and then two or three months of very high demand. For 
example, demand in the range of 50 to 100 would be received 
for several months and then a demand for 3300 would be re~- 
ceived. It would be ditficult to accept the assumption 
that the 3300 units demanded represented 3300 requisitions. 
The importance of this point is that if the 3300 units 
demanded were not on separate requisitions the value of 


unit effectiveness becomes distorted. 
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i PennOmmor RECORDING DEMAND 


Up to this point the problems associated with fore- 
casting demand have been defined as the ability to describe 
the probability density function of demand and the selec— 
tion of the proper forecasting scheme to predict the para-— 
meters of that density function. One of the major factors 
contributing to the complexity of these problems is the 
nature of demand patterns. An examination of the 1000 item 
sample used in this study indicated that in the majority of 
cases the demand patterns were very erratic. The primary 
effect of erratic demand 1S to cause large variances thereby 
rendering the task of fitting a probability distribution to 
PiemecnAang@mcata, particularily ditticuli ws More important, to 
the inventory manager, however, is the effect these large 
variances have on his inventory levels. 

One of the major factors cauSing demand patterns as the 
inventory manager’ sees them is the demand reporting system 
EhaitmwaSmutialized. fivamsship requires a repair part, then 
the demand is created at the shipboard level at a particular 
point in time. Under current procedures there are four 
echelons of supply from which the ship can obtain the repair 
part Bequired. Therefore*the initial point at which the 
demand 1S recorded depends on where the requirement 1S 
satisfied. The four echelons referred to are the ship 


itself, a mobile logistics support ship, a stock point and 
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finally, an inventory control point. These echelons are 
critical in determining when a demand gets recorded at 
the ICP, which is the inventory manager and has the res— 
ponsibility of setting system inventory levels. For 
example, if a shipboard demand occurs on the first of the 
month and the ship can satisfy the demand from its own 
stocks then the demand is recorded at the shipboard level 
only at that point in time. 

If the demand did not reduce the onboard stocks to the 
ship's reorder point this particular demand would not be 
recorded at the stock point until sufficient demands 
occurred to reduce the stock to the reorder point, at which 
time the ship would place a replenishment demand on the 
steck point. Another factor affecting the time at which 
the demand is reorded at the stock point is whether the 
ship 1S at sea or in port since the ship has no means of 
communicating (except for high priority requirements) 
directly with the stock point if it 1S at sea. The point 
here is that the demand that a stock point observes could 
represent demand for that day or an accumulation of demand 
over a period of time. Since stock points report daily 
to the ICP's this is also the situation that the ICP in- 
ventory manager faces. There is such a large population 
of ships that enter and leave port at different times it 
would seem that over a period of time, the system demand 
would smooth itself. However, a review of actual data 


patterns does not. substantiate this. Thus, an investigation 
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into the method of collecting data seems warranted. Since 
time would not permit a thorough investigation it was de- 
cided to look at the effects of this situation through a 
simple simulation model in which demands for a single item 
were generated from seven different ships. It was assumed 
that the time between demands is exponentially distributed. 
The time at sea waS assumed to be uniformly distributed 
between 0 and 45 days. Based on these assumptions the days 
on which demands occurred were randomly generated for each 
ship. Also the dates of arrival and departure from port 
were randomly generated, and initially all ships were con- 
sidered at sea. It was also assumed that all demands in- 
curred while the ship was at sea are reported to the stock 
point upon arrival in port. Additionally, any demands 
occurring while the ship was in port were reported to the 
stock point on the day they occurred. 

The alternative method of recording demand investigated 
in this simulation had the ship report the demand directly 
Bemthe LCP ron mshnesdayve1t O©Ccugrod.) Inherent in) this alter=— 
native is the assumption that the ship has the capability 
of reporting demand directly to the ICP which is recognized 
as unrealistic at this time. However, Since the ships 
currently keep records of daily demand through the 3M 
(Material, Management and Maintenance) system, it was felt 
that the above alternative was possible. Thus, we are com— 
paring consumption data vice demand data as currently 


reported in the Navy. 
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The simulation was run for a two year period with the 
demand recorded by month at both the ICP and stock point. 
Additionally, the mean and variance of demand over this 
two year period were computed for each reporting system. 


The results of this simulation are presented in Table VII. 


TABLE VII 


SIMULATED MONTHLY DEMAND 
AS RECORDED AT ICP AND NSC 


DATA YEAR MONTHLY DEMAND 


Consumption 
Data eee ee ey 7 eee 5 29.5 2 Bee 
Pole Slee ey ewe oe DOr 2 15 oS 
Demand 
Data fo eowmea 1) M6723 23° 20 “17 20. 


Poe 9. 13988 21582 .-4 40 08 22535 





Mean of Consumption Data = 15 Mean of Demand Data = 15 
Variance of Consumption Variance of Demand 
Data = 32 Data = 84 


As expected, the alternative of reporting demand direct— 
ly to the ICP (consumption data) shows a smoother pattern 
than that reported to the stock point (demand data). The 
critical result however is that the variance of demand at 


the stock point waS approximately 25 times greater than the 


variance of demand recorded at the ICP. With smoother or 
more regular demand patterns and smaller variances the task 


of forecasting future demands should be easier. It would 
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be presumptuous to draw any conclusions based on such 
tentative results. However, further investigation into 


the use of consumption data does seem appropriate. 


40 





V. CONCLUSIONS 


A. DEMAND FORECASTING TECHNIQUES 


The following recomendations are made with regard to 


the forecasting schemes to be utilized to predict the 


mean and variance of the lead-time demand. 


1) 


2) 


That the sample variance we be used for predic— 
ting variance of lead-time demand regardless of 
the probability distribution chosen. There 
appears no statistical justification for using 
the Mean Absolute Deviation regardless of the 
particular estimate of MAD used. Since there 

waS an apparent strong relationship between the 
sample MAD and the mean quarterly demand when 
ENCHEStUGy was comauctced) at FMSO am 1962 1t might 
be worthwhile to perform such a study with cur- 
rent demand. The use of the coefficients of 1.37 
and .717 computed in 1962 certainly does not 
appear warranted. 


The results of this study show no distinct pre- 
monCICCmmOL Lhemus: seid SanDlewdverdge Cimeeine 
exponentially smoothed estimate of the mean. 


Pee ROBABIIGITY DENSITY FUNCTION 


It is recommended that the current procedures of using 


three probability density functions to represent lead-time 


demand be replaced by the procedure of using the single 


distribution of gamma. This recommendation is based on the 


Following conclusions: 


1) 


Lace Gamnlda sex hioiteSPEnG@VverSaumlity Lo represent 

low demand items as well as high demand items. 
Additionally the gamma is a non—negative probability 
Gensity funetaon whoch icnsures thatvalll of the 
probability mass will be represented by actual 
demand (i.e. no such thing as a negative demand 

such as is the case with the normal). 
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2) 


$)) 


The gamma appears to be less sensitive to peaks in 
demand and therefore consistently assigns lower 
reorder levels creating much lower inventory hold- 
ing costs. The argument here is that unless the 
demand pattern iS smooth, stocking material under 
the normal assumption tends to put higher reorder 
levels than required to meet most demands. These 
higher reorder levels naturally result in higher 
protection levels but at the same time they require 
Eile syStem, C@ Camny a viot Of Steck and themefore 
incur shighn= holding #and investmemt costs £o Gretect 
against infrequent demands. 


Finally, the conversion from the present system to 
the gamma would require little effort Since the 
incomplete gamma function which gives the comple~ 
mentary cumulative distribution is already pro= 
grammed and tabulated. Additionally uSing the 
method of moments to compute the gamma parameters 
requires minimum machine effort. 
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APPENDIX B 


KOLMOGOROV-SMIRNOV TEST 
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* indicates gamma hypothesis is rejected 


44 





AP END he 
FORTRAN PROGRAMS ASSOCIATED WITH COMPUTATION OF INVENTORY 
COSTS AND UNIT EFFECTIVENESS 

FORTRAN programs used in calculating values for Tables 
mo andy. 

The three programs in this Appendix represent computa— 
tions uSing a variable Q. For the values on Table V the 
programs were modified to set Q = 1.0. Additionally these 
basic programs were run with a counter to compute the 
number of periods when shortages occur and the total number 
of periods when orders were placed. This later modifica-— 
tion was necessary to calculate the observed protection 
levels. 

Important labels used in the programs and not previously 


defined are as follows: 


BAR Mean monthly demand 
VAR Variance of monthly demand 
PT Desired protection level 


TSHT TOtd eins oS *ShOre 
THCS Total holding cost 


OUE Overammamiilome. feCcLiveness 
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FORTRAN PROGRAM FOR GAMMA (A,q) 
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APPENDIX E 
FORTRAN PROGRAM AND RESULTS OF SIMULATION MODEL 
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Total monthly demand recorded at Naval Supply 
56 


Total monthly demand recorded at Inventory 
Center (NSC) 


Day of departure from port for Ship I 
Control Point (ICP) 


Day Of arrivat in port for Ships! 
. Variance of demand recorded at ICP 
Variance of demand recorded at NSC 


Mean demand recorded at ICP 
Mean demand recorded at NSC 


Day of demand for Ship I 


Variable names used in the program are defined as follows: 
ISHIP (I) 

NPORT (1) 

NSEA (TI) 

NICPM 

NSCM 

NXBARI 

NVARI 

NXBARN 

NVARN 
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